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We present NLO QCD results for Wbb and Zbb production at the Tevatron including full bottom- 
quark mass effects. We study the impact of QCD corrections on both total cross-section and 
invariant mass distribution of the bottom-quark pair Including NLO QCD corrections greatly 
reduces the dependence of the tree-level cross-section on the renormalization and factorization 
scales. We also compare our calculation to a calculation that considers massless bottom quarks 
and find that the bottom-quark mass effects amount to about 8-10% of the total NLO QCD cross- 
section and can impact the shape of the bottom-quark pair invariant mass distribution, in particular 
in the low invariant mass region. 
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1. Introduction 

Precise theoretical predictions for the production of a weak gauge boson (W /Z) in association 
with a bb pair are very important both in the search for a Ught Standard Model (SM)-like Higgs 
boson and in the search for single-top production at hadron colliders, in particular at the Tevatron. 
Indeed, W/Z + bb represent the major irreducible backgrounds to the associated production of a 
light Higgs boson with a weak gauge boson (WH and ZH with H bb), which are the main search 
modes for a SM-like Higgs boson at the Tevatron |2|]. At the same time, Wbb is an irreducible 
background for single-top production, which is being measured for the first time at the Tevatron as 
pp ^ tb,tb with t{t)^b{b)W 

The signal cross sections for both WH/ZH and single-top production are known including 
higher order QCD and electroweak corrections. It is therefore important to also control the back- 
ground to a good level of theoretical precision. In the present experimental analyses the effects 
of Next-to-Leading (NLO) QCD corrections on the total cross-section and the dijet invariant mass 
distribution of the Wbb and Zbb background processes have been taken into account by using 
the MCFM package [Q], which implements the zero bottom-quark mass (m^ = 0) approxima- 
tion [| 0, |. 

In this proceedings we present results for the NLO QCD cross-sections and bb-pair invariant 
mass distributions at the Tevatron, including full bottom-quark mass effects [^. We separately 
study the impact of NLO QCD corrections and of non-zero bottom-quark mass effects. NLO QCD 
corrections stabilize the theoretical prediction of total cross-sections and distributions, reducing 
the dependence on the renormalization and factorization scales. On the other hand, the presence 
of a non-zero bottom-quark mass mainly affect the cross-section in the region where the bb-paiT 
invariant mass is small, both at Leading Order (LO) and at NLO in QCD, and amount to an overall 
8-10% difference with respect to the zero bottom-quark mass approximation. 

2. NLO calculation 

The hadronic production of a W boson with a bb pair occurs at tree level in QCD via the 
qq' Wbb partonic process. On the other hand, the hadronic production of a Z boson with a bb 
pair consists, at the tree level in QCD, of two partonic channels, namely qq — > Zbb and gg Zbb. 
The NLO QCD corrections to the tree level partonic cross-section consists of both one-loop 0{as) 
virtual corrections and C?(o;.v) real corrections corresponding to the emission of one extra parton 
from the tree level parton processes. The NLO hadronic cross-section is obtained by convoluting 
the parton-level NLO cross-sections with NLO Parton Distribution Functions (PDF). 

The 0{as) virtual corrections to the partonic cross-section contain ultraviolet (UV) and in- 
frared (IR) singularities. The UV singularities are calculated using dimensional regularization and 
cancelled by introducing a suitable set of counterterms (see Ref. for details). IR singularities are 
isolated using dimensional regularization and cancelled against the analogous singularities arising 
in the 0{as) real corrections to the partonic cross-section. 

The 0{as) virtual corrections to the partonic cross-section consist of one-loop self-energy, 
vertex, box and pentagon diagrams. We apply techniques developed in the NLO QCD calculation 
of Htt |[T^, |n|] for the calculation of scalar and tensor loop-integrals. In particular, we calculate 



2 



NLO QCD corrections to Wbb and Zbb production 



Laura Reina 



the tensor loop-integrals via Passarino-Veltman reduction (PV) [12]. We encounter instabilities 
due to the quasi-vanishing of the Gram determinant(s) of the process only in one box diagram and 
in several pentagon diagrams. We are able to obtain stable numerical results by combining, when 
necessary, sets of gauge invariant diagrams. 

The 0{as) real corrections to the partonic cross-section consist of the qq' Wbb + g and 
qg{qg) —>■ Wbb + q{q) subprocesses for Wbb production and of the qq — > Zbb + g, gg —>■ Zbb + 
g, and qg{qg) Zbb + q{q) subprocesses for Zbb production. We have extracted both soft and 



collinear IR singularities by implementing a Phase Space Slicing method with two cutoffs [ |13[ ] 
in order to isolate the soft (4) and collinear (5c) singularities respectively. Both the soft, hard- 
coUinear, and hard-non-collinear parts of the cross-section depend on the cutoffs, but their sum is 
cutoff independent over a large range of values of the cutoffs (see Ref. [Q]). 

Both analytical and numerical results for the NLO hadronic cross-section have been checked 
with two independent calculations based on different programming languages and public/in-house 
packages. The analytical reduction of the calculation has been obtained using FORM [Q and 
Maple codes, while the numerical results have been obtained using Fortran and C codes. The FF 
package [|l^] has been used to check some of the IR-finite scalar and tensor integrals. Finally, the 
hard non-colUnear real corrections have been double-checked using Madgraph [|I|,[l3, HI]. 



3. Numerical results 

In these proceedings we present results for Wbb and Zbb production at the Tevatron including 
NLO QCD corrections and a non-zero bottom-quark mass fixed at mfo = 4.62 GeV. Results for 
Wbb have been published in Ref. [^], while results for Zbb are currently being cross-checked and 
should be considered as preliminary. Both W and Z boson are considered on-shell and their masses 
are taken to be = 80.41 GeV for the Wbb runs and Mz = 91.1876 GeV for the Zbb runs, 
while, in each case, the other weak gauge boson mass is calculated via the relation Myv = cos Qw^z 
with sin^ = 0.223. The LO results use the 1-loop evolution of a, and the CTEQ6L set of 
PDF [|l|], while the NLO results use the 2-loop evolution of a, and the CTEQ6M set of PDF, with 
a^^^{Mz) =0.118. The W boson coupling to quarks is proportional to the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements. We take V„d = V^v = 0.975 and V„, = Vcd = 0.222, while we 
neglect the contribution of the third generation, since it is suppressed either by the initial state quark 
densities or by the corresponding CKM matrix elements. 



We implement the kj jet algorithm [20, 21, E3, 23] with a pseudo-cone size R = 0.7 and 



we recombine the parton momenta within a jet using the so called covaiiant ^-scheme []21[]. We 
checked that our implementation of the kj jet algorithm coincides with the one in MCFM. We 
require all events to have a bb jet pair in the final state, with a transverse momentum larger than 
15 GeV {p^j^ > 15 GeV) and a pseudorapidity that satisfies [t]'^'''| < 2. We impose the same pj 
and I T] I cuts also on the extra jet that may arise due to hard non-coUinear real emission of a parton, 
i.e. in the processes W /Zbb + g or W /Zbb + q{q). This hard non-coUinear extra parton is treated 
either inclusively or exclusively, following the definition of inclusive and exclusive as implemented 
in the MCFM code |j5|]. In the inclusive case we include both two- and three-jet events, while in 
the exclusive case we require exactly two jets in the event. Two-jet events consist of a bottom- 
quark jet pair that may also include a final-state light parton (gluon or quark) due to the applied 
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recombination procedure. Results in the massless bottom-quark approximation have been obtained 
using the MCFM code 
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Figure 1: Dependence of the LO (black solid band), NLO inclusive (blue dashed band), and NLO exclusive 
(red dotted band) total cross-sections on the renormalization/factorization scales, including full bottom-quark 
mass effects. The l.h.s. plot is for pp Wbb and the r.h.s. plot for pp — > Zbb . The bands are obtained by 
varying both jiR and jXp independently between /io/2 and 4/Zo (with /Xq = mi, +Mv/2 for V = W,Z in the 
pp Wbb and pp — > Zbb cases respectively). 

In Fig. |I] we illustrate the renormalization and factorization scale dependence of the LO and 
NLO total cross-sections, both in the inclusive and exclusive case. The bands are obtained by 
varying both and jXp independently between Ho/ 2 and 4/io (with /Iq = m/, +My/2 for V = W,Z 
in the pp — > Wbb and pp —>■ Zbb cases respectively), including full bottom-quark mass effects. 
We notice that the NLO cross-sections have a reduced scale dependence over most of the range of 
scales shown, and the exclusive NLO cross-section is more stable than the inclusive one especially 
at low scales. While the LO cross-section still has a 40% uncertainty due to scale dependence, this 
uncertainty is reduced at NLO to about 20% for the inclusive and to about 10% for the exclusive 
cross-section respectively. This is consistent with the fact that the inclusive NLO cross-section 
integrates over the entire phase space of the qg{qg) —>■ bbW /Z + q{q) channels that are evaluated 
with NLO as and NLO PDF, but are actually tree-level processes and retain therefore a strong scale 
dependence. In the exclusive case only the 2^3 collinear kinematic of these processes is retained, 
since 3-jets events are discarded, and this makes the overall renormalization and factorization scale 
dependence milder. This is better illustrated in Fig. ^for the case of pp Zbb, where the r.h.s. 
plots show the scale dependence of the cross-sections due to the single partonic channels (qq, gg 
and qg + qg). The strong residual scale dependence of the inclusive NLO cross-section is clearly 
driven by the qg + qg channel. Similar results are obtained for pp —>■ Wbb, as reported in Ref. [^. 

A first illustration of the impact of keeping a non-zero bottom-quark mass in the calculation 
of the NLO QCD cross-section is given in the l.h.s. plots of Fig. ^, where LO and NLO total 
cross-sections for pp Zbb are given, both for = and nih = 4.62 GeV, as functions of the 
renormalization and factorization scales (identified for the purpose of this plot). Neglecting bottom- 
quark mass effects overestimate the NLO cross-section by about 8-10%, depending on the choice 
of the scale. In Fig. ^ we analyze the impact of a non-zero bottom-quark mass on the bb-paii 



4 



NLO QCD corrections to Wbb and Zbb production 



Laura Reina 





5 




1 


— NLO massless _ 










MT O mticcivp 




4.5 






LO massless — 




4 






LO massive 










& 


3.5 








o" 


3 










2.5 


cuts; 


> 15 GeV 






2 




lTll<2 






1.5 




R = 0.7 

1 





0.5 



5 
4 

■B 3 

2 

1 






11/^0 



4.5 




1 


NLO massless ~ 

— NLO massive " 


4 






LO massless ^ 

LO massive 


3.5 








3 








2.5 


cuts: 


> 15 GeV 




2 




[r|l<2 
R = 0.7 




1.5 




1 




0.5 


1 


2 ^ 





1 ' 1 ' ' ' - 
— NLO massive 




qq initiated 




gg initiated 




qg initiated 




Inclusive case 








11/^0 



3 - 



S2 ^ 





1 ' 1 ' ' ' - 

— NLO massive 




qq initiated ~ 




gg initiated 




qg initiated — 




Exclusive case 






1 , 1 — i — i — i— 



0.5 



11/^0 



Figure 2: Dependence of the LO and NLO inclusive (upper plots) and exclusive (lower plots) total cross- 
section for pp Zbb on the renormalization/factorization scale, when ^Lr = jiF- The l.h.s. plots compare 
both LO and NLO total cross-sections for the case in which the bottom quark is treated as massless (MCFM) 
or massive (our calculation). The rh.s. plots show separately, for the massive case only, the scale dependence 
of the qq, gg and qg + qg contributions, as well as their sum. 



invariant-mass (m^j) distribution. We give results for both inclusive and exclusive distribution in 
the pp ^ Wbb case. In both cases distributions most of the impact is in the low m^g region, which 
can be important in a variety of different analyses. 
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